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A class of exact solutions of the Wheeler- De Witt equation for diagonal Bianchi type IX cosmolo- 
gies with cosmological constant is derived in the metric representation. This class consists of all the 
"topological solutions" which are associated with the Bianchi type IX reduction of the Chern-Simons 
functional in Ashtekar variables. The different solutions within the class arise from the topologically 
inequivalent choices of the integration contours in the transformation from the Ashtekar representa- 
tion to the metric representation. We show how the saddle-points of the reduced Chern-Simons func- 
tional generate a complete basis of such integration contours and the associated solutions. Among 
the solutions we identify two, which, semi-classically, satisfy the boundary conditions proposed by 
Vilenkin and by Hartle and Hawking, respectively. In the limit of vanishing cosmological constant 
our solutions reduce to a class found earlier in special fermion sectors of supersymmetric Bianchi 
Q\ • type IX models. 
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I. INTRODUCTION 



Three decades after the first exploratory steps were made |l],^| canonical quantum gravity continues as a vigorous 
program of fundamental research. New life has been breathed into this field about a decade ago by Ashtekar's 
^ ■ discovery [EJ of a representation of general relativity in terms of new variables, which render the Hamiltonian and 
t^- ' diffeomorphism constraints more tractable. The developments in nonperturbative canonical quantum gravity which 
C^l followed this advance justify a reasonable hope that a mathematically consistent quantum gravity might be 

attainable. Still, the program is far from completed and many facets of the theory remain to be explored. One 
such facet is the relation between the metric representation and Ashtekar's representation. This question is far from 
trivial: Ashtekar's formulation starts out with complexified general relativity and suitable reality conditions on the 
new variables have to be imposed at the end. The metric representation, on the other hand, stays within the domain 
of the real theory all along. Thus it is not clear, a priori, whether there is a unique one to one relation between the 
complexified theory and the real theory. Can, e.g., a single solution of the Wheeler-DeWitt equation in Ashtekar's 
variables before imposing reality conditions give rise to several mathematically and physically distinct solutions in 
metric variables? 

In the present paper we wish to examine this question in the framework of the minisuperspace model of diagonal 
Bianchi type IX with a non- vanishing cosmological constant. 

Kodama || and Brencowe ||lfj[| have found a simple solution of the basic constraints of quantum gravity with 
cosmological term in Ashtekar's variables in the form of an exponential of the Chern-Simons functional. By projecting 
C3 ' to Bianchi type IX geometries one also obtains a solution for the minisuperspace model || . 

In the present paper we shall start from the Wheeler-DeWitt equation for the minisuperspace model in the metric 
representation and a specific choice of operator ordering (which is gleaned from the special operator ordering appearing 
naturally if supergravity is used as a starting point [jll] [l4j]). The Ashtekar representation is then introduced only 
on the quantum level as a mathematical device, like a Laplace-transform, to simplify the equations. In fact, the 
representation is introduced as a kind of complexified momentum-representation, in which Ashtekar's variables are 
the complexified canonically conjugate momenta of the inverse triad corresponding to the Bianchi type IX 3-metric. 
The integration contour in the complex manifold spanned by these momentum variables may be chosen quite freely 
within the requirements of convergence and the vanishing of boundary terms in partial integrations. Integration 
contours which can be deformed into each other while satisfying these requirements are topologically equivalent. 
However, a given solution in Ashtekar variables may admit topologically inequivalent choices of integration contours. 
Such a solution in Ashtekar variables may then correspond to several mathematically and physically distinct solutions 
in the metric representation. In fact, we shall show that this happens for the Chern-Simons topological solution in 
the diagonal Bianchi type IX minisuperspace model. We find that five topologically inequivalent integration contours 
over the Ashtekar variables exist, which are organized by five distinct saddle-points of the reduced Chern-Simons 
functional and the accompanying paths of steepest ascent and descent. These findings raise the interesting question, 
whether similar results may also be obtained in the full theory. The answer to this question is not obvious, because 
the enlargement of the configuration space, in principle, could render integration contours topologically equivalent, 
which appear as topologically inequivalent when projected on the minisuperspace under investigation. 

While this general question transcends our minisuperspace framework and must be left open here, our results for the 
minisuperspace model yield several new exact solutions in metric variables with non- vanishing cosmological constant. 
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These solutions turn out to be of interest in their own right. To simplify their discussion we restrict ourselves to 
the physically more interesting case A > throughout this paper and postpone the examination of the case A < 
to a future work. We discuss the asymptotic limits of the solutions for h — > and A — > and show that, at least 
semi-classically, two of them satisfy the cosmological boundary conditions proposed by Vilenkin jl5| and by Hartle 
and Hawking [ |l6| , |l7| |, respectively. Furthermore, we show that it is just the no-boundary state which additionally 
fulfills a physically well-motivated normalizability condition. 

The remainder of this paper is organized as follows: In section |n| we establish our notation, put down the Wheeler- 
DeWitt equation in the adopted operator ordering, and briefly list the five exact solutions for vanishing cosmological 



constant which are known from recent work on Bianchi type IX supergravity. In section III the Ashtekar repre 



sentation is introduced as a complexified momentum-representation, and the resulting Chcrn-Simons solution of the 



Wheeler-DeWitt operator is given. In section IV that solution is transformed to the metric representation along five 
topologically distinct integration contours, which establishes a basis of five linearly independent solutions which are 
all generated by the Chern-Simons solution. Their asymptotic behavior for fi — > 0, A — > 0, or A a 2 — > oo (where a 
is the scale parameter and A the cosmological constant) is studied in section [v| Here, in addition, we establish the 
relation of two of these solutions to those picked out by the boundary or no-boundary conditions of Vilenkin, and 
Hartle and Hawking, respectively. In the last section we draw some conclusions and indicate how our results may also 
be used to establish certain limiting forms of Bianchi type IX models coupled to a scalar matter field with very small 
mass. 

II. METRIC REPRESENTATION 

A. Wheeler-DeWitt equation 

The purpose of this section is to establish some notation and to derive the Wheeler-DeWitt equation for the Bianchi 
type IX model in the metric representation in a specific factor-ordering. We start from the Einstein-Hilbert action 
with a cosmological constant A 

S EH [g^] = J d 4 x^j (k-2A) (2.1) 

M 

where a possible boundary term has been omitted, since such a term will not contribute to the resulting Lagrangian 



density. In (2.1) the action integral is taken over the 4-dimensional space-time manifold of the Universe, while g and 
1Z are the determinant and the Ricci-scalar of the 4- metric g= (g^ u ), respectively. Performing the ADM space-time 
split a ^apse-function N, a shift-vector N l and the 3-metric h%j of the spatial slice are introduced in the usual way 



split a tapse-tunction i v , a stujt- vector iv 
fl8| . Then the action (|3]) takes the form 



At / d 



j3. 



X ■ 



1G 



— ( 3 ft - K 2 + KijK ij - 2A) , (2.2) 



where h = det {hij), 3 1Z is the curvature scalar of the spatial manifold, and Kij is the extrinsic curvature tensor. The 
spatial homogeneity manifests itself in the existence of infinitesimal coordinate transformations x H — x % + e J £(j)( x )> 
which leave tensors on the 3-manifold form-invariant jl9| . In the Bianchi type IX case the algebra of the Killing- vectors 
with (x) = g|j is chosen to be 

£(*)> = £ ^ k £( fc ) ' 

implying do) 1 = ^ Sijk & Aa} fc for the invariant basis Cj 1 [ p9| . Using this special basis in the following, the components 
of all tensors on the 3-manifold become functions of time t only. Let us now consider the Bianchi type IX form [T^-|2^1 

Mt) = e 2 " (t) (e mt % (2-3) 

for the 3-metric, but with restriction to the diagonal case. Then (/Jy) is a diagonal, traceless 3 x 3-matrix, which can 
be parameterized as 

C%) = diag(/3+ + V3 - V3 /3_,-2/3+) . (2.4) 



The action (2^2) yields a Lagrangian 
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C = 2 nNe 3a ( 3 U- K 2 + K i3 K ij - 2A) , (2.5) 

where the spatial integration has been carried out using Jii'A^Aii 3 = 2 (47r) 2 , and 3 1Z and ify have to be expressed 
in terms of the new metric variables {a, /3±} and their time-derivatives. It turns out, that C is independent of N and 
N l , so the conjugate momenta satisfy the primary constraints 

n » := ?£=0 n i := ^L = 0} (2.6) 

ON dNi 

i.e. they vanish identically. The preservation of the primary constraints in time leads to the secondary constraints 

^ = 0=^ (27) 
dN dW { <) 

the first of which is the Hamiltonian constraint, while the second one, in the present case, is solved by taking N % = 0. 
Therefore the 4-metric is now of the form 

g = -N 2 (t)dt 2 + h i:j (t) . (2.8) 

A straightforward Legendre transform yields the Hamiltonian 

N p~ 3a iii 
H= — H , where H = —n a + tt+ + tT + e ia U(B±) + 3 (8tt) 2 A e 6a . (2.9) 

48 7T 



Here the momenta 



ir a = ^ , tt ± = ^S- (2.10) 

da d(3 ± 



have been introduced, which are connected to the generalized velocities as follows: 

Furthermore, the curvature-potential U(/3±) has been defined as 

U(P±) ■= -6 (4tt) 2 e 2a 3 K = 3 (4tt) 2 Tr (e 4/3 - 2 e~ 2/3 ) . (2.12) 

The Hamiltonian constraint now simply reads H = 0. 

To quantize the model, one may seek for solutions of the Schrodinger equation 7i |W) = E where Tt is now 
interpreted as a self-adjoint operator over a suitably defined Hilbcrt space of wavefunctions |^). The Hamiltonian 
constraint, which must be satisfied by the physical states, then implies a restriction to zero-energy states, yielding the 
Wheeler-DeWitt equation 

H\V)=0. (2.13) 

The quantized version of the Hamiltonian ( |2.9| ), which is obtained from the usual rules of canonical quantization, 
suffers from the well-known ambiguity in the choice of th e factor-ordering. We shall resolve this ambiguity in such a 
way that a class of simple semi-classical solutions of ( 2.13| ) (to be given in the following section) , whose existence and 



form, due to their semi-classical nature, are independent of the factor-ordering, become exact solutions. Thi s w ill be 



achieved through the following non-standard procedure: It is easily checked that the classical Hamiltonian (2.9) can 
be written in the formQ 



1 These definitions, and the resulting factor-ordering, are suggested by the existence of a supersymmetric extension of this 
form of the Hamiltonian [ pd|Jl^ , |l^ | . The factor-ordering chosen here is not contained in the class considered in p4| . 
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H = 



in a - $ a 



Mr" + #, a 



in + — $ . 



in + + $ 



in — $ 



in + <E> 



+ 3 (8n) 2 A e 6a , (2.14) 



where $ is defined as 



$ := 27re 2a Tre 2/3 



(2.15) 



Assuming now canonical commutation relations [a,7r Q ] = [(3±,n ± ] = ih, the momenta in the {a, /3±}-representation 
may be expressed as n a = —ih d a , n ± = —ih d± , and H becomes 

H = \ hd a - $ a 1 [ hd a + <f>, a ]-\ ?i<9 + - $ J [ + $J - [ hd- - $ _] [ ft<9_ + $ _] + 3 (8tt) 2 A e 6Q . (2.16) 



Finally, this corresponds to a Wheeler-DeWitt equation 

j?i 2 [9 2 -a 2 -d 2 _ ] - 12h$ + e 4a U(0±) + 3(8n) 2 Ae 6a }^\n. i :A) n 



(2.17) 



in the metric representation. The non-standard term ~12?i$ is a quantum correction to the classical "potential" 
e 4a U(j3±) and appears as a result of our choice of operator ordering in equation ( [2.16] ). 



B. Solutions without cosmological constant 



For completeness, and because they will play some role in the discussion of the general sol ution s in the case A ^ 0, we 
shall now present a derivation of the five known solutions of the Wheeler-DeWitt equation ( 2.17 ) without cosmological 
constant fn^Q. 

An obvious, exact solution of the Wheeler-DeWitt equation with vanishing A can be extracted from the form ( 2.16| ) 
immediately [0|: 



<™ := cxp 



(2.18) 



In the classical limit h — > one can interpret S = i$ as the Euclidean action of this wavefunction. As known from 
Hamilton- Jacobi theory the derivatives of the action with respect to the generalized coordinates play the role of the 
generalized momenta, so the classical trajectories may be computed via (2.11) from 



da 
iNdt 



-3a 



24tt 



d/3 



± 



3 -3a 



i Ndt 24tt 



(2.19) 



Clearly, because $eR, no real solutions to ( p. 19 ) exist, but one can introduce a new parame ter of imaginary time 



dr := iNdt and look for Euclidean solutions, corresponding to positive definite 4-manifolds via (2.8Q. It can be shown 
that the solutions of ( p. 19] ) in this Euclidean regime form a two-parameter family of classical Universes, which beco me 
all asymptotically flat and isotropic in the limit a — > oo, where a :— 2e a is the average scale factor. In fact, ( 2.18) ) is 
the well-known "wormhole-state" of the Bianchi type IX model. 



To derive other solutions of the model without cosmological constant we will first subject the Hamiltonian (2.16) 
to the similarity-transformation 



H = e 



(2.20) 



which yields 



H' = h 2 



dl~d\-d 2 _ 



2h 



*, Q d a - * + d+ - $ _ cL 



3 7r 2 a 6 A 



(2.21) 



Moreover, new variables playing the role of the inverse triad are introduced 



a i: --a 2 e- ft (>0) , ie{l,2,3}, 



(2.22) 
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where fa are the diagonal elements of the matrix (faj). Denoting derivatives with respect to these new variables with 
di := Tj^- and rescaling the cosmological constant into A = one finds for the Hamiltonian in the di-representation: 



H' = 12H 2 H 'i : where H 'i : = VjVk 



djd k - di + ^ a 



(2.23) 



Two important features of H' in this representation should be mentioned: 

• H' is a sum of three terms, producing each other by cyclic permutation of the indices. 

• H' is invariant under converting two of the <Ji into their negatives. 

The first of these properties suggests to look for solutions by solving the reduced equation 



0102 



d\d 2 - d 3 



= 



(2.24) 



with a ^-symmetric function W'. The simplest solution to this ansatz turns out to be 
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exp (T\ + cr 2 + 03 



(2.25) 



which by multiplication with the wormhole-state gives a second solution, the "Hartle-Haw king state" Jl 2[— 11 4[| . This 
name can be justified by discussing the equations for the classical trajectories analogous to ( [2.19 ). There it turns out 
that each member of the two-parameter family of solutions describes a regular (and isotropic) Universe in the limit 
a — > 0, i.e. in this sense ^° NB satisfies the "no-boundary" proposal. 

Three further solutions to the Wheeler-DeWitt equation occur just because of the second property of the <7j- 
representation of H', namely jl4| 



:= 



exp <7i — <jj 



Sijk 



(2.26) 



These "asymmetric solutions" create classical Universes which turn out to be asymptotically three dimensional: In 
the limit a — *■ oo the i-th dimension is curved to zero, while the remaining manifold becomes an asymtotically flat, 
spatially two dimensional wormhole. 

We shall add a few remarks on these states concerning their normalizability. An investigation of the five wave- 
functions reveals that ^® VH and the three asymmetric solutions ^ are bounded functions on minisuperspace, whereas 
^5r B is unbounded (in fact, it grows super-exponentially for a — > +oo, when (3± is kept fixed and small). In the limit 
a — ► — oo, i.e. considering vanishing scale factors, all five solutions approach unity, and thus a normalization integral 
over the full {a, /3±}-space diverges in any case. However, the four bounded solutions ^ VH and ^ may at least 
be called normalizable in the distributional sense. Introducing a suitable integration weight they will even become 
normalizable in the usual sense. 



III. ASHTEKAR REPRESENTATION 



A. Wheeler-DeWitt equation in Ashtekar's variables 



The transformation of the Wheeler-DeWitt equation ( 2.17 ) into the Ashtek ar re presentation |25|,[26j has been well- 
prepared during the last section by introducing the new triad variables <Ji via (2.22). All that remains to be done now 
is to transform (2.23) to the generalized momenta A% conjugate to <jj by performing a suitable Fourier-transformation. 

It will be crucial for all that follows that we shall choose not a standard Fourier-transformation which is carried 
out along the real axes, but a generalized, complexified Fourier-transformation defined asQ 



¥(a) = J d 3 A e 3 ' X *'(A) 



(3.1) 



2 Since we are dealing with complex- valued quantities A% anyway, the standard "i" has been absorbed in these new variables. 
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Here we assume that £ C C 3 is a smooth, three dimensional manifold with the following properties: 



(i) The integrand of (3.1) and its first and second derivatives with respect to Ai vanish at the border <9E of E. 

(ii) E is completely contained within a domain where ^'(A) is a holomorphic function of Ai- 

Assuming these properties of E it is possible to convert ai and di to the new variables: It turns out that they obey 
exactly the classical conversion rules 



di -» Ai 



-di 



(3.2) 



where for brevity di := jrj-, i.e. the form of the Fourier-transformed equation is independent of the choice of E and 



unaffected by our complcxific ation procedure! 
The transformed equation (2.23) now reads 



H'& = 0, where H' = £ H[ , H[ := djd k A 3 A k - Ai - - di 

i=l L 

which we will refer to as the Wheeler-DeWitt equation in the Ashtekar representation. 



(3.3) 



B. Solution with cosmological constant 

An exact solution of the Wheeler-DeWitt equation with non-vanishing cosmological con stant can now be easily 
constructed: As in section II. B we try to make vanish each term of the Hamiltonian of (3^3) seperately, i.e. we seek 
for a solution obeying the three equations 



Hi & = d,jd k 



2 



§' = 



(3.4) 



Furthermore, one may even try to solve 



= o 



(3.5) 



for each set of {i,j, k} — {1, 2, 3} simultaneously. This requirement gives the unique solution 



= exp 



~ ( 2A!A 2 A 3 -A 2 



(3.6) 



This solution is even known for the general spatially inhomogeneous case, where the exponent is given by the Chern- 
Simons functional [p||To|] . In pa the semi-classical content of the wave function ([[!]) was analyzed within the Ashtekar 
representation. Here we shall be interested, instead, in its transformation to the metric representation. In j^] the 
transformation back to metric variables was also attempted, but without success. I n th e following, we shall make 
use of the freedom in the choice of contours in the generalized Fourie r-tra nsformation (3.1) to derive several different 
solutions in the metric representation which are all generated from ( |3.6| ) by the choice of topologically inequivalent 
contours. 



IV. TRANSFORMATION TO THE METRIC REPRESENTATION 



A. General form of the transformation 



By def mtion , the generalized Fourier-transform (3T) of ( |3.6| ) represents an exact solution of the Wheeler-DeWitt 
operator (2.23) in the Oi-representation: Each wavefunction of the form 



^'(«r) = / d°A exp 

s 



I V' R ~ \a 2 + AiA 2 A z 



(4.1) 



where S is chosen to imply a sufficient ly l arge fall-off for the integrand on t?£, is a solution for the Bianchi type IX 
model with cosmological constant. In (4.1) new variables 



Ki := — A Gi — — A a 2 e @* 
2 12 



(4.2) 



have been introduced. We will show that there exist several manifolds S satisfying the conditions (i) and (n), corre- 
sponding to different solutions in the metric representation, so obviously, due to the existence of several topologically 
inequivalent contours, the gene ralized Fourier-transformation is not unique. It will turn out that for the space of 
exact solutions defined via (4.1) the number and location of the saddle-points of the integrand's exponent 



F(A, k) := A- iZ-^A 2 + A1A2A3 



(4.3) 



will play an essential role. While for the asymptotic form of the solutions obtained from ( }4.1| ) by the saddle-point 
method the importance of the saddle-points is obvious and well-known |27j , p8fl , their importance for the exact solutions 
is a surprise, which arises because of the freedom in the choice of the integration contours. These saddle-points are 
determined by the equations 



dF 
dAi 



— = & m - Ai + AjAk = 



Sijk 



1 



(4.4) 



which in the case K\ ^ k 2 A3 ^ ±1 may be rewritten in the form 

A\ = «i + -42-4,3 , 

A 2 (1 - .A3) = k 2 + ki A 3 , 
(A 3 - k 3 ) (1 - A 2 3 ) 2 = (Kl + K 2 A 3 )(/t2 + Kl A3) ■ 



(4.5) 
(4.6) 
(4.7) 



ird equatio n ( 
points via ([4.(^ 



4.7) is of fifth order and so yields five different solutions, which correspond to five different 
and (O). 



The th: 
saddle- 
One may show that also in the remaining case K\ = k 2 (4.4) has five different solutions with A 3 -components still 
obeying (4.7). The fact that there always exist five saddle-points of the integrand's exponent will be shown to result 
in a five dimensional space of solutions ( [4.lj ). 

For concreteness, let us assume a special representation of sucessive integrations 



= / dA 3 / dA 2 / dAi exp 

r 3 r 2 (^ 3 ) r 1 (A2,A 3 ) 



A 



F(A, k) 



dA 2 exp 



- ( 2A 2 k 2 + 2A3K3 - A 2 2 - A 3 



r 2 (-4 3 ) 



dAi exp 



-( -Af + 2Ai (K1+A2A3. 



ri(^ 2 ,^ 3 ) 



(4.8) 



which one may show to be of no restriction to the general case. The one dimensional Gaussian Ai-integral just has a 
single saddle-point located at 



Ax = Ki + A2A3 



(4.9) 



Up to a factor that may depend on A, but that will be absorbed in a proportionality sign "cx" in the following, there is 
only one non-trivial value this integral can take: Since the integrand has to vanish at the ends of Y\ each integration 
curve can be deformed into the curve of steepest descent, which in the new coordinate 
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A[ =Ai- (ki + A 2 A 3 ) 



(4.10) 



turns out to be simply the real axis. So the Ai -integration yields 

+00 



J dA[ exp ~ [~A^ + («i + A2A3 



\tt exp 



j (k 1 +A 2 A 3 ) 2 



(4.11) 



and (4.8) is turned into 



*'((?) oc J dA 3 exp 



k\ + 2„4 3 k 3 - A\ 



r 3 



&A2 exp 



A 



(1-^) + 2A 2 (k 2 + ki^ 3 



r 2 (-4 3 ) 



The ^-integral is again of Gaussian form with just one saddle-point 

K 2 + Kl-4 3 



l-Al 



(4.12) 



(4.13) 



To keep this saddle-point away from infinity, one has to exclude the points A3 = ±1 from the integration path r 3 . 
With this prescription and using the new coordinate 



X = A 2 \ 1 - A 



K 2 + K1A3 



the integration is again easily carried out along the real axis with the following result: 



d-4 3 



VT—A. 



exp 



- ( k{ + 2„4 3 k 3 - Ai 



(k 2 + K1A3) 2 
l-Al 



(4.14) 



d-4 3 



exp 



A V 1 - A 3 I + A3 



2 A 3 K3-AI 



(4.15) 



Here further variables n± := (ki ± n 2 ) have been introduced. Let us now define the new exponent as 



/(*,*) 



1 - z 1 + z 



2 z K3 — z , 



(4.16) 



then ty' is easily expressed as 



®'{a) cx J 



dz 



VT^z* 



exp 



(4.17) 



Furthermore, a new coordinate u = Arcsin z will prove useful, yielding the representation 

1 



ty'(cr) (X J du exp 



A 



/(sin u, k) 



(4.18) 



where now the points u = ±-| (and all 27r-periodic repetitions) have to be exclu ded from the new integration path C. 
It is possible to show that all saddle-points of / are determined by the equation (4.7) when A3 is replaced by z and the 
z-so lutions are translated to the u-plane afterwa rds. S o the five saddle-points of the original integral representation 
(4.1) occur in the one dimensional representation (4. IS), too, and we have still all freedom to choose a specific solution. 
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For a discussion of the location of the saddle-points we shall now restrict ourselves to the case A > 0. Then it turns 
out that three of them are always real, which we will denote by 



< -1 , -1 < z + < , z 3 > 1 



(4.19) 



in the following. The other two solutions may also be real or conjugate complex depending on the values of the real 
positive variables Kj. The two corresponding regions in the K-space (which is the minisuperspace with fixed A) are 
separated by a caustic which is characterized by the existence of a marginal saddle-point of f(z,K). We will refer to 
the first part of the minisuperspace, where all five saddle-points appear real in the z-plane, as the "Euclidean regime" 
and denote the two additional real saddle-points byn 



< 



< z < 1 



(4.20) 



whereas in the "Lorentzian regime" there exist two complex saddle-points labelled according to the signature of their 
imaginary parts as 



3z V ii < > 



(4.21) 



The caustic defined above will play an important role for some particular solutions which get their dominant integral 
contribution at the corresponding marginal saddle-point; for fig. |l| it has been computed numerically in the {k, /3±}- 
space (here n = Aa 2 , see below) . If on the other hand solutions are considered which receive contributions from 
several distinct saddle-points, it should be clear that the caustic is of no significance to them. However, such solutions 
will turn out to be of little physical interest anyway. 




FIG. 1. The caustic in minisuperspace for A > 



Knowing the saddle-points of /, one may now calculate the curves of steepest descent for the two different regimes 
by solving the equations 



:/(sinu(r) ) = 3/(sinw e ) , ge{-,+,vu,0,3} 



(4.22) 



The result of a numerical approach to this problem is given in fig.g. In the stripe |5Ru| < ^ all solutions of (4.22) are 
presented. The dashed curves reach to +oo with respect to / and are given just for completeness. It is remarkable 
that there exist paths running into the singularities at ^3 = ±1, corresponding to z = ±1 and u = ±5, in such 
a manner that 5R/ tends to —00. This results in additional possibilities to create integration contours obeying the 
requirements of the generalized Fourier-transformation. 



The naming of the indices introduced here will be justified later in the discussion of the asymptotic behaviour of the solutions. 
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Euclidean regime 



Lorentzian regime 





Saddle -point * Singularity at u = ± 



SR/ decreases 



FIG. 2. Saddle-points and curves of steepest descent 



B. Basis of linearly independent solutions 



With the knowledge of the curves of steepest descent now a basis of linearly independent solutions (4.1) may be 
defined by choosing the following integration paths: 




FIG. 3. Basis set of integration paths 



The same paths can be chosen in the Eucl idean and in the Lorentzian domain. It is easily seen that all curves end at 
essential singularities of the integrand of ( 4.1 8| ) in such directions that the integrand and all its z-derivatives vanish. 
This guarantees that 



^' e {(y) cx J du exp 



/(sin u, k) 



ge{-, +, vu , 0, 3} 



(4.23) 



are indeed solutions of the Wheeler-DeWitt equation in the a-representation. Furthermore, any path satisfying the 
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fall-off condition can be deformed into a superposition of the curves defined in figure ||, consequently one may express 
each wavefunction in terms of the corresponding solutions (4.23). Thus the above mentioned basis property of ( 4.23| ) 
is proven. 

For concreteness, let us write 



* e (d0 = V° WH Af e (X) / du exp 



A 



/(sin u, k) 



(4.24) 



where Af g (\) is a normalization factor that usually is defined via (\I r e |\f r £ 



1. Here 



denotes the scalar product 



of the underlying Hilbert space which we do not know explicitly. But even if we knew it, an analytical treatment 
of this normalization procedure would certainly not be tractable. That is why we adopt another convention and 
normalize the wavcfunctions by the condition 



* e (a = 0) 



(4.25) 



which implies that "J e becomes independent of A for a — > 0. This requirement gives analytically soluble normalization 
integrals, yielding 



JV+ = A/1 = Af = JV 3 = - 



2ie tl a^o 2ie~ 



2e^ 



A^O 1 



■ e >■ 



2TrI (n)+iK Q {Li) 



Att 



with ti := — 
' 2A 



(4.26) 



where K and Jo are the McDonald's and the modified Bessel function with index 0, respectively. The asymptotic 
behavior of A/" e in the limit A — > will be useful in the next section. 

There are other solutions of the Whceler-DcWitt equatio n wh ich will turn out to be rather interesting, and w hich 
of course can be expressed in terms of the basis solutions (4.24): By choosing another order of integrations in (4.8) 
the wavef unctions "J 1,2 defined by 



#i(ki,k 2 ,k 3 ) == *3(k2,K3,«i) and *a(«l)«2,«3) :=*3(«3»«i 5 «2) 



(4.27) 



apparently solve the Wheeler-DeWitt equation, and the asymptotic behavior in section ^ will show that the relations 

*iH • *->h ,; + ' Hill (4-28) 



hold in the limit Aa — > 00. Since (4.24) is a basis with non-vanishing ^> g in this limit, it follows that (4.28) must 



hold in general. This claim is proven by considering a general expansion in terms of the basis states, with coefficients 
depending on A, but not on a, where one makes use of the limit a — > 00 at fixed A =/= 0. 
With the solutions \&i and ^2, now the sum 



^3 = \I>i + ^2 + *3 



(4.29) 



turns out to be a solution as well, which is symmetric with respect to arbitrary permutations of the Kj . If in addition 
one uses the fact that the manifolds 



£+ := 



I Ae<D 3 I AjeRe M %, je{l,2,3} } 



(4.30) 



which are symmetric under permutations of the Aj , are suitable to perform the generalized Fourier-transformation in 
the representation (4.1), sucessive integrations analogous to (4.8)-(4.15) reveal that ^0 and are symmetric under 
permutations of the Kj . 

Finally, by examining the integrand's symmetries under complex conjugation, ^ e turns out to be a real wavefunction 
for q e {-, +, 0, 1, 2, 3} , whereas VEV;; is complex-valued in general. 
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V. ASYMPTOTIC RESULTS 



In this section we will discuss the asymptotic behavior of the solutions to the Wheeler-DeWitt equation w ith n on- 
vanishing cosmological constant der ived above. This will be done by evaluating the integral representation (4.24) in 
the limits h — > and A — > 0, where ( 4.24 ) turns out to assume the form of a saddle-point integral. Consequently, the 
asymptotically leading term of such an integral 



dz 



.*/(») 



A^O 




2vrA 



.*/(«.) 



(5.1) 



will be of particular interest |27| , Pq| . Here z s denotes that saddle-point of / which provides the only contribution to 
the integral in the limit A — > 0, and the sign of the square root has to be adjusted to the direction in which this 
saddle-point is passed through. In the following u ^/~" will denote the principal value of the square-root. 
The application of (5.1) to the general representation (4.24) yields 



exp 



x /(sinu) 



2 du' 



■/(sinu) 



Afg(X) V^A 



exp 



■/(*) 



(/"(Z)(Z2_1) + ^( Z ) Z ) 



ge{-,+,v«,3} 



(5.2) 



Since A = M, this formula includes both, the limit ?i — > and A — > 0. That the curves defined in fig. ||| give indeed 
the saddle-point contributions mentioned in (E>.2|) can be extracted from fig. ^ by analyzing the topological properties 
of / in detail. The asymptotic expansion ( |5.2| ) also holds for ^> , as long as just the Euclidean regime is considered. 
In the Lorentzian case all saddle-points can be passed through by Co and may th erefo re contribute, so in dependence 
on the variabl es a and /3± one has to choose the highest saddle-point to employ ( |5.2| ). 
The result (5^2) now suggests the definitions 



S„ 



A e := 



\^"{z e ){zl-l) + f{z e )z Q 



A^O 



/ 7rA J\f g A e exp 



n 



(5.3) 



In the limit fi — > the exponent S can then be interpreted directly as the action of the wavefunction (up to a 
constant term, which may arise from A/" e ), while A plays the role of Gaussian fluctuations around the saddle-point z B . 
Considering the limit A — > 0, these interpretations hold no longer; here S and A may just be called the phase- and 
amplitude function of 'J, re spec tively. 

Unfortunately, equation (4.7) determining the saddle-points is of fifth order, so analytical expre ssion s for the roots 
are not available. This is why we shall first restrict ourselves to the isotropic case (3± = 0, where (4/7) can be solved 
explicitly, yielding 



zf = -1 



4 0) = i 



(o) 



1 



1 



1 



Jo)* 



K < 



K > 



where K := — Aa 2 . (5.4) 



With these preparations we shall now turn to a detailed discussion of the specific solutions. 



A. Vilenkin state 



1. Semi-classical limit Ti — > 



If the saddle-point Zy}, is inserted into (5.2) one obtains the action 
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o(0) 



3 in 
~A~ 



1-1 



a 2 A 



(5.5) 



for the isotropic case. Choosing the lapse function N = 1, the only non-trivial equation for the classical trajectories 
reads 



da 
~dl 



1 dS^ 



(0) 



3a7r 9a 



1 - -a 2 A 
3 



(5.6) 



and is easily integrated to give 

a(t) = ao cosh(aQ 1 t) , where ao 



A ' t>0 



(5.7) 



Reinserting this classical DeSitter solution into (5.6) and remembering that 11 denotes the principal value of the 
square-root reveals that, as mentioned, (5.6) is solved only with the restriction to t > 0, i.e. collapsing Universes are 
not described by ^ vu . 

Metrics with a < ao can be obtained by solving the Euclidean version of ([5.6]) with dr = i di, yielding 



a{T) = a sin (a rj , — - < r < a 7r 



(5. 



The restriction of the r- variable appears for the same reas on a s discussed above for (5.7). Denoting the line element 
of the unit 3-sphere by dO, the 4-metric corresponding to (5J3) reads 



ds 2 = a 2 , ( dr' 2 + sin 2 r' dtt 2 ) , with r' := 1 r . 



(5.9) 



It describes exactly a 4-half-sphere with radius ao , which in the limit r' — ► ^ may be extended to the DeSitter solution 
( |5.7p . Furthermore, the point a = 0<^r' = 7risa regular point of the manifold, i.e. \t vjl satisfies the no-boundary 
proposal in the isotropic case. 

It is interesting to see whether these properties remain true when anisotropic corrections are considered. The 
calculation of such corrections is straightforward: One has to expand z vil for small (3±, and insert this expansion in 
the expression for S vil . One finally obtains^] 



S vil = Svu + 367ri a 2 



1 - § a 2 A 



24 + a 2 A 



{fil+fil)+0{(3l) . 



This action yields for the classical trajectories in imaginary time in leading order of (3± (with n — Aa 2 ) 



da 

d7 



-V 1 -4k + 4 



(« + 2) 2 



(k + 3) 2 - 10 
V 1 - 4k 



(5.10) 



(5.11) 



d/3± _ 2 3 + V 1 -4k 
dr a 2 + k 



/3± , 



(5.12) 



where we have assu med k < j. Since the prefactor of /3± in (5.12) is positive definite, the anisotropy decreases with 
decreasing t, while ( 5.11 ) then tells us that a tends to qq. Consequently, a flat, cylindrical 4-geometry 



ds z = dr 2 + a 2 dn 2 > 



(5.13) 



is approached in this limit with ao playing the role of the cylinder radius. 



4 A corresponding expression for Svu, which is here obtained as a limit of the exact result, was first derived by Del Campo 
and Vilenkin, using the WKB method |^]. These authors did not discuss the consequences for the semi-classical trajectories 
generated by this action. 
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With increasing r the point a = is approached and the anisotropy grows rapidly so that the validity of (5.11) and 
(5.12) breaks down. Therefore a = is not a regular point of the Euclidean space-time manifold any more, and the 
no-boundary proposal is not fulfilled for ^> vil . 



A discussion of a Lorentzian version of (5.11) and (5.12) reveals that classical DeSitter-like Universes are described, 
which grow exponentially in time t while the anisotropy decreases monotonously.^] However, in general j3± does not 
tend to zero, but approaches a finite value. 

In any case \EV;i describes an expanding Universe, i.e. quantum mechanically speaking ^ vu supports a current in 
minisuperspace which is directed to the positive a-axis. Thus this wavefunction satisfies the condition proposed by 
Vilenkin and so is identified as the Vilenkin state of the Bianchi type IX model with cosmological constant. 



2. The limit A -> 



In contrast to the limit h — > 0, in the ca se A — > the location of the saddle-points themselves depends on A via the 
variables Kj. Nevertheless, the expansion ( |5.2| ) remains applicable as long as just the leading term is considered, and 
happily now the A-corrections of the saddle-points can be calculated taking account of the full influence of anisotropy. 
The first terms read 



Zvu = \ <73 A + - (rio-2 A 2 + C(A 3 ) 



(5.14) 



and the phase- and amplitude functions are calculated to be 



S va = i$-i?if (of+of + of) +0(A 2 ) , 



A vu =1 + 4- (of + a! + of) + 0(A 3 ) 



(5.15) 



Using (5.3) and the asymptotic behavior of J\f g in accordance with (4.26), one finds for the wavefunction 



A^O 



Af Vil V^X A vu e5 Svu 



A^O 



(5.16) 



Thus in the limit A — > the wormhole state of the A = -model is approached. 



3. The limit k — > 00 



Writing the wavefunctions ^ g defined in (4.24) in the alternative form 



* e = K-hK / du ex P 



YI-kk 2 ( sin u + cosh [2-y/3/3_] 



hA 



-2/9+ _ I ^ Sinu 



smu 



,2/3 4 



(5.17) 



another possibility to approximate a saddle-point integral occurs, namely the limit k = Aa 2 — > 00. This includes 
the cases A — > 00 and, in particular, a — > 00. 

An asymptotic expansion of the saddle-points is again possible without any restriction for the anisotropy variables, 
but the result is rather lengthy. So we turn at once to the expressions for S vil and A vil obtained in this limit: 



Svu = ^ J - 4 v^ 3 + ^ Tr (e- 2 ? \ e 



4/3 



7 - Tr e 2/3 • Tr e 



Tre 



6/3 



(5.18) 



i 1 , 
77 ^Tre 2 
2 



+ 0(k _ ») 



(5.19) 



5 In the Lorentzian regime the action calculated from (5.3) has a non- vanishing real- and imaginary part for f3± 7^ 0, so here 
the definition of classical trajectories is not a priori clear. Since we are interested in pseudo-Riemannian 4-geometries, and 
because $tS V u dominates QSvu for large scale parameters a anyway, we choose the real part of S V u to discuss the Lorentzian 
classical trajectories. 
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Surprisingl y, the contribution i$, which usually arises from the wormhole state and would be expected to give a term 
of O(k) in ( |5.18 ) has completely disappeared. 

To shorten the final expression for \& v «, let us now expand for small mean anisotropy 8 := J ' 2 + (3 2 : 



Vn 



^° K o (&) -2^(11) 



12/3 2 



-|f(6/3) 2 
e exp 




ITT 

T 



(5.20) 



From this result it is clear that |\I/ v , ii | 2 is bounded for a — > oo. Furthermore, a saddle-point expansion for 8 — > oo 
at fixed a reveals that \&vs! is square-integrable over B±, and because it is bounded for a — ► — oo (cf. (4.25)) it is 
normalizable in the distributional sense.^J 

As a general result the Vilenkin state becomes concentrated about 8± = in the limit k — > oo, but with a 
non-vanishing Gaussian width 



(5.21) 



To give an idea of the behavior of the exact analytical solution we have computed the real- and imaginary part of the 
wavefunction numerically in dependence on a and A, assuming 8± = and picking units with fi = 2n Q 



9*% 





mm 




FIG. 4. The Vilenkin state 



The conspicious oscillations, which start on the caustic k = i Aa 2 = 3, indicate that the initial Euclidean action 
has turned Lorentzian, corresponding to a pseudo-Riemannian classica l Un iverse. For A — > the Gauss function of 
the wormhole state is recovered, whereas our normalization condition (|4.25| ) is responsible for the form of the graph 
at a — 0. 

In addition, the value of \^ vu \ 2 C on the caustic has been computed for fig.|5|in dependence on the anisotropy variables 
0± at fixed values A = 3, h = 2ir. Since all the Lorentzian trajectories end at the caustic, we follow Hawking and 
suggest to interpret |^| \ as the distribution of the initial values for a classical evolution of the Universe. For the 
Vilenkin wavefunction the distribution is nicely concentrated about 8± — 0, but there are directions of the 8±— plane 
in which l^valc takes a finite value for 8 — > oo. Moreover, one may show that in these directions a tube with 



This entails the possibility to use Marolf 's method |29| for introducing a scalar product to define a Hilbert space of physical 
states, which contains ^ V u, presumably as a ground state. 

7 This unusual choice of Ti is due to the use of differently scaled variables in the numerical work and is of course of little 
significance to the figure. 



15 



finite height and width is approached. This is not in contradiction to the normalizability of the wavefunction in 
the distributional sense, since on the caustic a — ► — oo in the same limit. However, it implies that |^ Vi! |g is not 
square-integrable with respect to the p± -variables! 




FIG. 5. The initial value distribution generated by the Vilenkin state 



In subsection C we will define a new state as a superposition of ^ vii and ^* vil in such a way, that this undesirable 
property of the initial value distribution is removed. 



B. Asymmetric states 

1. Semi-classical limit Ti —* 



In the isotropic case (3± = the actions of th e w avefunctions ty± and '5 3 are easily calculated as the corresponding 
saddle-point contributions in accordance with (5.4): 



S 



(0) c-(O) 



6ni 



1 + 5 



2 ao 



a 
2a 



(5.22) 



Obviously these actions are purely Euclidean for generic a and A, and this statement remains true if anisotropic 
corrections are considered. An additional discussion of the limit k — > 00 will show that the asymmetric solutions 
themselves are indeed of no physical relevance, nevertheless, an examination of the case A — > will be worthwile. 



2. The limit A -> 



As in the Vilenkin case, in the limit A — * asymptotic expressions for the saddle-points z± and Z3 are available 
taking full account of anisotropy. The results read 



z± = -l± - \a x -<r 2 | +0{\ 2 ) 



(5.23) 



and using them, the calculation of the phase- and amplitude functions is straightforward, yielding 



S+ = i 



— + * + ?i ( cr 3 =F I CTi — cr 2 1 ) 



■0(A) 



#3 = i 



67T 

— + $ - h ( cr 3 - ax - a- 2 



■0(A) 



(5.24) 



and 



A± = 2 



1 - 4 ( °"3 T \0\ - 0"2 I 



+ G(A 2 ), A 3 = - 



1+4(0-3- O-i - (72 ) 



-C(A 2 



(5.25) 



1G 



Thus one obtains for the wavefunctions in the limit A — > 



lim * 



± 



exp 



-a 3 ± \a\ - o- 2 1 



lim * 3 = f w exp 



A— 



CT3 — (7l — CT 2 



(5.26) 



where, in addition, the asymptotic behavior (4.26) of the normalization factors has been taken into account. Obviously 
t he so lutions ty± are not differentiable at a\ — a 2 , be. in particular at [3- — 0, but if they are replaced by \ti,2 via 
(4.28), solutions with nice analytical prope rties are obtained. These are related to ^3 by permutations of the Oj (or 
equivalently the mentioned in (KL27J). 

The wavefunctions ze{l,2,3} are now easily recognized as asymmetric states, which approximate the corre- 
sponding solutions ^ of the A = -model when A tends to zero. 



3. The limit k — > 00 



Let us first consider the case K\ > K2, where the saddle-points can be expanded as follows: 



K2 . K 3 [k\ - K\) 



K 2 K 2 K\ Kry 

The asymptotic behavior of the solutions ^± and ^ 3 then turns out to be 



+ e>(0 , z 3 = k 3 + 



Kl K 2 



h > 



+ 0{k,- 1 ). (5.27) 



\Ef u 

WH 



1 



K2,l K3 



«1,2 



1.2 



exp 



1 1 



^2,1 K 3 
«1,2 



(5.28) 



and 



K. — *■ OG 2 



^o(M) ^ K 3 



K X K 2 



exp 



A V k 3 



(5.29) 



Since f(z,K) is invariant under the permutation «! <-> k 2 , the correspo ndin g resul ts for the case K\ < k 2 can now 
easily be obtai ned b y exchanging the Kj-indices 1 «-> 2 in the equati ons ( 5.27 )-( |5.29| ). The asymmetric states and 
^2 defined in (4.27) then have to be constructed from ty± by using (4.28). Consequently, the asymptotic behavior of 
^i, i e {1, 2, 3}, may be written in the closed form 



2 VttX 



KjK-k 



exp 



.2 r» K j K k 



where £j 7 -fc = 1 



(5.30) 



We should mention that these solutions diverge badly in the limit n — > 00, namely like e Aa , so they are surely not 
normalizable in minisuperspace for any sensible choice of the scalar product. That is why we reject them as candidates 
for the physical quantum state of the Universe, and we are left with a just two dimensional, physical space of solutions 
spanned by 9 vil and ^y iV 



C. The no-boundary state 



In this section we will show that there exists a superposition of \t vi! and ^* vill such that this wavefunction is 
normalizable in minisuperspace and square-integrable on the caustic. Moreover, this uniquely determined solution 
will turn out to satisfy the no-boundary condition proposed by Hartle and Hawking Jl6| , pTt (at least in the sense 
that the classical Universes described by this wavefunction are regular at a — 0) f\ One may construct this solution 
by normalizing the Vilenkin state to approach unity in the limit f3+ — * +00 at (3- — on the caustic. If one then 
considers the difference of this new Vilenkin solution and its conjugate complex solution, the obtained distribution 



A qualitative discussion of a no-boundary state for the anisotropic Bianchi IX metrics was first given in |20| and, more 
explicitely, including a numerical plot for the wavefunction in the semi-classical limit, in An expansion for small anisotropy 
has already been given in Ji(J. However, the results obtained there are very lengthy and hard to interpret. For a discussion of 
the semi-classical trajectories generated by the no-boundary state in Ashtekar's variables see also J2(| 
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on the caustic is obviously square-integrable with respect to (3± (for further explanation of this construction cf. fig. 
Finally, we shall choose the still unspecified overall normalization factor as usual in accordance to our convention 
.25[). The solution defined by this procedure turns out to be 



* ™ := *° 



° WH 3 I Af NB J du exp 



(5.31) 



where 



M NB = 2e/ 



J {v)-iK {v) 



K (iJL)Jo(v)-2nIo{ti)Ko(v) 



1 4 
' M= 2A' vt= \ 



(5.32) 



Here the special integrals 



2 OO 

Jo(y) '.= J Ax exp [— v sin 4 x\ , fCo(v) := J Ax exp \—v cosh 4 x\ 



(5.33) 



have been introduced which, to a certain extent, may be considered as generalized modified Bessel functions.^ 

It is clear from its construction that ^ NB is integrable on the caustic, normalizable in minisuperspace in the distri- 
butional sense (as ^ va ) and we shall show below that it satisfies the no-boundary condition for h — ► 0. Furthermore, 
^ NB obviously is a real-valued wavefunction. The beha vior o f the no-boundary state for k — > oo can be immediately 
extracted from the asymptotics of the Vilenkin state ( 5.20|) . In this way an asymptotic description of $ KB in the 
Lorentzian regime is available, so we will restrict ourselves to the Euclidean regime throughout the following. To 
discuss the limits Ti — > and A — * we shall first expand the normalization factor M NB for small A. Saddle-point 
expansions of J7q and /Co m the corresponding limit v — > oo finally yield for the leading term in the asymptotic series 



a — o 2e» a^o 
N NB ~ — == - iN , 



(5.34) 



where we have also used ( 1.26| ) . Obviously the normalization factors J\f NB and Afo have the same asymptotic behavior 
in the limit A — ► 0, where ^ NB may be written in the form 



du exp 



A 



/(sinu, k) 



(5.35) 



If one now chooses the integration path C_ © C+ © C3 Co, which is equivalent to C vil © C* u , the final expression for 
the asymptotic behavior in the limit A — > becomes: 



,0 1 




(5.36) 



This representation displays nicely the individual saddle-point contributions in the Euclidean regime and will prove 
useful for the following discussions. 



1. Semi-classical limit Ti — > 



Since for h — > the saddle-point zq always provides the dominating contribution in comparison with z± and Z3, the 
relation (5.36) implies 



9 Similar normalization integrals occur in the calculation of Mvu and Ao in ( |4 ■ 2€)| ) , but with squared trigonometric functions 
in the exponent. Such integrals lead to the modified Bessel functions / an d Kg , that alternatively may be expressed as 



hypergeometric functions of the i_Fi-type. As a generalization, the integrals ( |5.33| ) may be written in terms of generalized 
hypergeometric functions of the ^F^-type, but the integral /Co requires logarithmic contributions (as Ko) which, as far as we 
know, have no special name in the 2f2-case. That is why we prefer to deal with the integral representations themselves. 
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(5.37) 



and all that remains to be conside red is an expansion of in the limit h — ► 0. 

Using once more the expansion ( |5.2| ), the action of the solution is calculated to be 



S, 



(0) 



3ni 




(5.38) 



i n th e isotropic case, where also the normalization factor Ao has been t ake n into account (in contrast to defintion 
( |5.3|) ). The expression is easily recognized as the negative Vilenkin action ([Tf|). Consequently the classical trajectories 
and spacetime metrics are the same as in the Vilenkin case up to a reversal of the r-direction and, as there, the no- 
boundary condition is satisfied in the isotropic case. 

But let us now consider the influence of anisotropy: Then the action is of the form 



36ni a 



1 



24 + a 2 A 



{Pl+pl) + 0{(3l) 



(5.39) 



implying (with k = Aa 2 ) 



dr 



4k - 



(k + 2f 



(k + 3) 2 



10 



V 1-4k 



(5.40) 



d/3± 
dr 



2 3 - V 1 - 4k 



0± 



(5.41) 



for the Euclidean, classical trajectories. As the prefactor of f3± in (5.41) is po sitive definite, the point @± = is 
attractive for decreasing r, i.e. in this r-direction j3 tends to zero. Then in (5.4C) the p±-term may be neglected, and 
the scale factor will reach a = at a finite value of r, say r = 0. The asymptotic form of (5.40) and (5.41) for r — > 
simply reads 



do 
d7 



1 



d/3± r^o 2 
dr a 



fti 



(5.42) 



and the j3± -equation can be integrated to give 



dr 



In &1 



>o 2 da 
a dr 



dr 



lna^ 



^ (3± oc 



(5.43) 



Consequently, the classical Universes become exactly isotropic in reaching a — 0, and there look the same as in the 
case (3± = 0. Therefore also Universes anisotropic at a > remain regular as a tends to zero. So \fr NB is indeed a 
solution of the Wheeler-DeWitt equation which satisfies the no-boundary proposal semi-classically and its name is 
justif ied, a fter alL_ 



If ( 5.40 ) and ( 5.41 ) are considered with increasing r, the anisotropy grows exponentially as in the Vilenkin case, and 
the validity of these equations breaks down. Finally, the trajectories approach the caustic, and tunneling processes 
to Lorentzian space-time become more and more probable. 



2. The limit A -> 



To discuss the limit A — > of the no-boundary state ^ NB we will again make us e of the relation (5.36). Since 
— > for A — > is known from before for ie{l,2,3}, only the behavior of in ( 5.36| ) remains to be discussed. 
An expansion of for small A is straightforward and we just give the final results: 



(5.44) 



So 



— + $ - n(ai + (7 2 + CT 3 ) 



+ G(A) 



(5.45) 
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A 



1 + - (<Ti + (T 2 + 03) 



C(A 2 ) , 



(5.46) 



lira * 

A^O 



1 T 



exp 



(71 + 0"2 + (73 



= $ u 



(5.47) 



Obviously the no-b ounda ry st ate (2.25) of the A = -model is approached by \&o- From equation (5.36) it now 
follows, using ( |5~26l ) and ( fj| ) 



lim \P, 

A^O 



5 (g^ -*-) 



(5.48) 



We can see that $ KB does not approach the no-boundary state of the A = -model when A tends to zero 
even though both wavefunctions (as many others!) satisfy the no-boundary condition semi-classically for h — > 0. 
The difference comes from our additional requirements that ^S NB for A ^ be normalizable in the distributional 
sense with respect to all three variables {a, (3±) and be square-integrable with respect to (3± on the caustic^} The 
normalizability of ^ NB for a — ► 00, like that of ^ vu , is directly related to the fact that, semi-classically, these states 
describe a Lorentzian Universe in that limit. E.g. had we defined the no-boundary state simply by v&o this would have 
given an acceptable limit \I>° B for A — > 0, but for A 7^ the semi-classical limit of "Jo gives, besides a Lorentzian, also 
Euclidean contributions. As a consequence, this alternative "no-boundary state" would contain components, which 
describe additional Euclidean Universes for a — > 00, and make the wavefunction diverge in this limit. Furthermore, it 
is easily checked that "to is a wavefunction which is not square-integrable on the caustic. 




A 6 



FIG. 6. The no-boundary state 

Figure ^ shows a numerical plot of the no-boundary state in the {a, A}-plane analogous to fig.|| To allow a clear view 
the plot has been bounded to values of lying in the interval [—5, +5]. Like in the Vilenkin case, rapid oscillations 
start in the Lorentzian regime, and clearly ^> NB — ► 1 as a — > 0. In view of the strong increase of the amplitude 
with a at small A, one might worry about normalizability in minisuperspace. However, in connection with (5.31) the 
asymptotic result for the Vilenkin state Q5.20 ) guarantees that 5™ falls off as for each A if sufficiently large scale 
factors are considered. 



10 The property of normalizability in minisuperspace opens again the possibility to define a scalar product for physical states 
as in P9|, and to consider the no-boundary state as a ground state in the resulting Hilbert space. 
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FIG. 7. The initial value distribution generated by the no-boundary state 



Finally, fig. [7] shows the resulting initial value distribution |\&jvb|c for the classical evolution of the Universe. As 
announced above, I'FnbIc rapidly approaches zero with increasing j3 in all directions of the /3±-plane, thus representing 
an integrable probability distribution. The figure shows that the quantum state ^ NB implies an isotropic classical 
Universe in the Lorentzian regime, apart from possible quantum fluctuations about /3± = 0. 



VI. CONCLUSION 

The central purpose of this paper was to derive exact solutions of the Wheeler-DeWitt equation with cosmological 
constant by employing a generalized Fourier-transformation to the Chern-Simons solution (3.6). As a result we found 
five linearly independent states, which approach the five known A = -states when A tends to zero. However, for 
A ^ 0, there is just a two dimensional subspace which appears to be physically relevant in the sense that the states are 
normalizable wavefunctions in minisuperspace and allow a classical Universe with pseudo-Riemannian geometry. We 
were able to construct a basis of this subspace in form of the solutions ^ vit and ^ NB , which, semi-classically, satisfy 
the boundary conditions proposed by Vilenkin and Hartle and Hawking, respectively. Moreover, with the additional 
requirement of integrability on the caustic, just one solution remains, namely the no-boundary state $ KB . We may 
add a few remarks comparing the two physically interesting solutions $ m and ^ vil . 

If the Euclidean geometries arc considered, which arise from the semi-classical trajectories in the Euclidean domain 
of minisuperspace by introducing an imaginary time variable dr = i dt, the two solutions behave very differently: 
While a Universe described by the no-boundary state becomes isotropic and remains regular in approaching a = 0, 
the Vilenkin-trajectories diverge in some (3± -directions for a — > 0, i.e. the corresponding Euclidean 4-geometry behaves 
singular for small scale parameters. 

If one considers the Euclidean approach to the caustic then all trajectories generated by the Vilenkin state reach 
the caustic in the isotropic point k = 4, (3± = 0, whereas the no-boundary trajectories run to infinite anisotropy. 

The creation of a classical Universe with pseudo-Riemannian geometry happens by a tunneling process through this 
caustic. The tunneling probability is given as the value of on the caustic and thus is a function of f3± only (A taken 
as fixed). The results for these distributions, which rule the initial values for the classical evolution of the Universe, 
are shown in fig.|^ and fig.0 and predict both an isotropic Universe as the most probable classical Universe among 
all alternatives. However, just the no-boundary distribution is sharply concentrated about (3± = and integrable 
with respect to the anisotropy variables, whereas the Vilenkin distribution is broader and not integrable (even though 
bounded). 

In the Lorentzian regime both wavefunctions describe a DeSitter-like Universe, which becomes more and more 
isotropic as a tends to infinity. As a difference, always describes expanding and collapsing Universes simultane- 
ously, while ^ vu just describes the expanding phase, and thus is not invariant under time reversal. 

In summary one might prefer the no-boundary solution ^ NB as the quantum state of the Universe. However, we 
should stress the fact that it is not the no-boundary idea itself which picks out uniquely this solution among the five 
solutions described by the Chern-Simons functional, but several additional, physically well-motivated integrability 
conditions have to be imposed on the solution to get $ KB . Then, the resulting wavefunction is found to satisfy the 
no-boundary proposal in the semi-classical limit h — > 0. 



Finally, it may be of interest to consider the application of the results obtained here to a model with a massive, 
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scalar matter-field)^]. We want to show in brief that the A-solutions found above are directly applicable to this model 
in the limit m — > 0: 

A massive Klein-Gordon field <f> gives a contribution 



S m at [g^, tf>] = J d^xy/^g (g^ '<j> ^<j) >v + m 2 (j) 2 ^ 



M 



(6.1) 



to the action fl2.1|) and finally yields an additive term 



*2 q2 i 2 6 2 j.2 

-a ca + 7T a m 



(6.2) 



in the Hamiltonian (2.16). If now the field amplitude <f> is rescaled with the mass m, rj :=m4>, H mat is converted to 



*2 2q2 , 2 6 2 

-n m o' + irarrf 



(6.3) 



Comparing with (2.21) it is easily seen that in the limit m — > this term corresponds exactly to the contribution 
which arises from a cosmological constant A = ?/ 2 . Thus one obtains the asymtotic solution for the Bianchi type IX 
model with a massive, scalar field 



lib — ' u "? 9 \ 

*(a,/3±,^;m) ~ \&(a,/3±; A = m 4> ) 



(6.4) 



at fixed A = m 2 2 , which, by further expansion about m = 0, may be extended to small, but non- vanishing masses 
of the field. 

As an interesting project for future consideration we leave the application of the generalized Fourier-transformation 
to the general form of the Chern-Simons functional |^,|l0|. It would be interesting to determine whether again a 
Vilenkin and a no-boundary state are obtained as topologically inequivalent functional integrals from the Chern- 
Simons solutions. Work in this direction is in progress. 
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